X-ray absorption branching ratio in actinides: LDA+DMFT approach 
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To investigate the x-ray absorption (XAS) branching ratio from the core 4d to valence 5/ states, 
we set up a theoretical framework by using a combination of density functional theory in the local 
density approximation and Dynamical Mean Field Theory (LDA+DMFT), and apply it to several 
actinides. The results of the LDA+DMFT reduces to the band limit for itinerant systems and to the 
atomic limit for localized / electrons, meaning a spectrum of 5/ itinerancy can be investigated. Our 
results provides a consistent and unified view of the XAS branching ratio for all elemental actinides, 
and is in good overall agreement with experiments. 

PACS numbers: 78.70.Dm, 71.10..W, 71.27.+a 



Understanding the physics of elemental actinide sohds 
is an important issue for many body physics as well as 
for apphcations in nuclear power generation. In the early 
actinides, the / electrons behave as waves delocalized 
through the crystal, while in the late actinides, the / 
electrons behave as particles localized around each atom. 
Plutonium is near the localization-delocalization edge 
separating these two regimes. X-ray absortion (XAS) 
from the core 4c? to the valence 5/ in conjunction with 
atomic physics calculations has been a powerful probe 
of the evolution of the valence and the strength of the 
spin-orbit coupling across the actinide series. 

A large number of spectroscopies have been applied to 
this problem. For example, photoemission spectroscopy 
of Pu, has revealed a multiple-peak structure in the oc- 
cupied part of the density of stateaii^i^. A combination 
approach of local density approximation (LDA) and Dy- 
namical Mean Field Theory (LDA+DMFT) has allowed 
the interpretation of these features in terms of / electrons 
which arc delocalized at low frequencies with a mixed va- 
lent electron count^. Other interpretations of Pu spectro- 
scopies using LDA+DMFT with other impurity solvers 
have been presented recentlji^i^ii. 

High-energy probes such as electron energy-loss spec- 
troscopy (EELS) or x-ray absorption spectroscopy (XAS) 
constitute a different set of spectroscopies which have 
been intensively used to study the electronic structure 
of the actinide series^i^iiiSiiiii^i^ii^. These experimental 
works, combined with theoretical calculations exploiting 
a powerful sum rule, and the electronic structure of the 
atom has yielded valuable insights on the degree of local- 
ization of valence 5/ electrons in actinide o^'^'^° as well as 
the spin-orbit strength^ i^°'^^ . 

Here we address the computation of the branching ra- 
tio within DMFT. There are several motivations for this 
study. 1) While the atomic multiplet approach of Ref . [l^ 
describes very successfully the majority of the data in 
the late actinides, it is restricted to the case that the / 
electrons are strictly localized. It is therefore useful to 
extend this approach by embedding it in a more general 
method that captures the itinerant limit as well. 2) Many 
other spectroscopies of the actinides such as photoemis- 
sion spectroscopies are not well described by a localized 
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FIG. 1: Schematic diagram of the XAS from core 4d to va- 



lence 5/ transition, 
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tions, respectively. The core levels are clearly discretized by 
4^3/2 and 4^5/2 level, but the valence levels are overlaped 
between 6/5/2 s,nd 6/7/2 level due to mltiplet splittings and 
electron itinerancy. 



model. LDA+DMFT provides a unified picture which 
reconciles the results of high-energy and low-energy spec- 
troscopies. In physical terms, DMFT can describe the 
low-energy part of the excitations as itinerant, and the 
high-energy excitations as localized. Different spectro- 
scopies weight different parts of the spectra. It is useful 
to incorporate the XAS branching ratio in the general 
LDA+DMFT formalism to achieve a unified interpreta- 
tion of high energy and low energy spectroscopies. 3) 
The branching ratio can be expressed as a ratio of two 
quantities, one involving the spin-orbit coupling and the 
other involving the / occupation. Therefore a proper 
theoretical interpretation of the experiments requires the 
evaluation of these two quantities. Until now, experi- 
mentally the occupancies were determined to be in an 
interval of allowed values. XAS and EELS experiments 
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were used to constrain the variations of the spin-orbit 
coupling for that range of occupancies. 

Figure [T] shows a schematic diagram of the Ad 
5/ transition in XAS. K. Moore and collaborators 
have shown that EELS experiments give equivalent 
informatio n^^'^^ . One envisions a large splitting due 
to spin-orbit coupling of the core levels. The electric- 
dipolc selection rule (Aj = 0, ±1), gives rise to two 
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(4(i3/2 — > 5/5/2). The branching ratio B is defined by 

S = 15/2/ (4/2+ /3/2). 

When the spin-orbit splitting in the valence band is 
negligible in comparison with the different sources of 
broadening (due to multiplet splittings or electron itin- 
eracy) shown schematically in Fig. [1] we can neglect the 
spin-orbit splitting in the final states and the branching 
ratio B is given by a statistical value of 3/5, which re- 
flects the relative degeneracies of the initial states. Hence 
the branching ratio is a probe of the strength of the spin- 
orbit coupling^ii^iii. 

Detail calculations based on the isolated atomic model 
were shown to be consistent with the experimental result 
of the late actinidcs. However, deviations get larger when 
the system becomes delocalized as in U and Np. Due to 
the delocalization, the spin-orbit strength is apparently 
suppressed. We will show that although LDA is a proper 
approximation for describing many properties of itiner- 
ant system, it underestimates the experimental results 
of early actinides. So more realistic theoretical model, 
which can treat both the atomic and itinerant physics, is 
needed. The LDA-I-DMFT method is the most promising 
method for this purpose. 

The spin-orbit sum-rule has been fully derived in 
the pioneering work of van der Laan, Thole and 
collaboratorsi2ii2ii2iiiSj and applied to numerous systems, 
(for recent review, see Ref. llSl ) Here we summarize the 
main steps of their derivation, indicating the places where 
the atoms are considered in a solid state environment, 
and where a treatment going beyond band theory and 
atomic multiplet physics, such as Dynamical Mean Field 
Theory (DMFT) is required. 

The X-ray absorption resulting from a core-valence 
transition, is described by a term in the Hamiltonian that 
couples the electromagnetic field to a transition operator 
Tq. which in the electric dipole approximation is given 
by: 



(1) 



Tg{i) is the atomic operator at site i. We consider 
transition from a core level denoted by jc (= Ic^s) to the 
partially occupied valence levels denoted by jy (= l-u ± s) 
with absorption of q polarized light. 



^ {lys;jymy\rq\lcs;jcmc)f^j^^^dj^,n^, (2) 



where dj^m^ is the annihilation operator of a core electron 
and fj^rn„ is the creation operator of a valence electron. 

The matrix element can be given as a product of a 
reduced matrix element containing a radial integral and 
an angular dependent part^^. 

{lys; jymy\rq\lcs; jcrric) = (-1)''""^'= [^t,^cjc]^^^ 

^'^ ^ -^^^^ ^ t ){l.s\\r\\hs), (3) 
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where [a, 6, • ■ • ] is shorthand for (2a-|-l) (26-1-1) • ■ • . Here, 
we used the Wigner-Eckart theorem. Because we are 
interested in the transition of given core (Ic) and valence 
(ly) states, we will omit the reduce matrix element below. 

When we consider the isotropic spectrum, where the 
light polarizations are averaged, the absorption intensity 
summed over the final states |/) from a many-electron 
ground state \g) is 

~E E E ^'^\a'm'J^^™-'vf]^.m^'^3am^) 

X {lcS\ jcm'^\Vq\lyS] jym'^){lyS\ jymy\Vq\lcS-, icrric) ■ (4) 

Here we assume that the interference terms {{Tq{i)Tq{j)), 
i ^ j) are negligible which is a valid assumption when 
the core electron states are effectively localized. Since 
\g) does not contain holes in the core level, the core shell 
operators are removed by d| ^, \g) = 0, and the creation- 
annihilation term in the intensity can be obtained as 



:) ~ (/iuini /J„m„ )'^™c™e (^) 



Where we neglected the core valence hybridization and 
we neglect the core-hole valence interaction, which is a 
reasonable approximation in Ad — > 5/ transition.— 

Combining Eqs. (3)-(5), the the dipole transition prob- 
ability is 



xy( 1 ) ( -'^ 1 
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=Mc.c]{i: I (6) 

We have used the normalization condition. 

Jc 1 jv \ ( jc 1 j'y 
^\mcq -my I \ nic q ~m'^ 

rncQ 

= [jv]^^Sj^_j'^Sm^m'^ (7) 



We consider the dipole transition from core 4c? states 
{jc = 2 ± 1/2) to valence 5/ states (j„ = 3 ± 1/2). Using 
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the relevant values of 6-j symbol in Eq. (5), the absorp- 
tion intensity for each transition is given as 

4/2 - /(4d3/2 ^ 5/5/2) = (ng/2)(2Z. + - (8) 



4/2 = ^(4rf5/2 ^ 5/5/2) + /(4d5/2 ^ 5/7/2) 

= K/2)/4 + (4/2)(24-l). (9) 
The branching ratio for Ad — > 5/ transition is given as 



B = 
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("7/2) + K/2)/[4(2/. ^ 1)] 
("7/2) + ("5/2) 



(10) 



For the valence states of orbital angular moment I, the 
expectation value of the angular part of the spin-orbit 
interaction is related to the electron occupation numbers 
(71 ) of the total angular momentum levels j± = / ± 1/2. 



]=j± 



j{j + 1) - l{l + 1) 



= -{l + l){n,_)+l{n,^). (11) 



Using the definitions (71'') = (715/2) + ("-7/2) Eq. 
pip , the spin-orbit sum rule for Ad 5/ transition is 
given by 



15 (77.'') 



(12) 



While we have made the assumption that the core elec- 
trons are localized, this assumption is not necessary for 
the valence f electrons. Therefore under the conditions 
stated before, the sum rule is valid not only in atomic sys- 
tem but also in solid system, and the spin-orbit strength 
can be estimated from the partial occupancy of valence 
states at a given site, theoretically. It also shows that the 
angular part of the spin-orbit strength is linearly related 
to the branching ratio in XAS and EELS^^'^^J^. 

To obtain the spin-orbit strength and branching ra- 
tio, we calculated the partial occupancies (rij^) in solid 
system by using the LDA and LDA-f-DMFT method?^ 
We use a relativistic version of the linearized muffin- 
tin orbital (LMTO) method for LDA calculations^^. 
In LDA-t-DMFT method, the itinerant spd electrons 
are treated using LDA, and the strongly correlated 
/ electrons are considered in DMFT approach, which 
maps a lattice problem to a single impurity problem 
in a self-consistent electronic bath^ii^^. To solve the 
impurity problem, we used the vertex corrected one- 
crossing approximationSl^ and the results are further 
cross-checked by the continuous time quantum Monte 
Carlo method^!. The Slater integrals F''{k = 2, 4, 6) and 
spin-orbit coupling constants are computed by Cowan' 
atomic Hartree-Fock (HE) program with relativistic 
corrections^^. We scale the Slater integrals by 70% to 
account for the screening of the solid. We take Coulomb 
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FIG. 2: The LDA expectation value of the angular part of 
the spin-orbit interaction. The LDA results are denoted by 
circles for actinide elements as a function of the 5/ electron 
count (n/). Paramagnetic phase is considered in all cases. 
Corresponding experimental data are denoted by the same 
color of dashes. The three common angular momentum cou- 
pling schemes are shown: LS, jj, and intermediate coupling 
schemes. 



interaction U = 4.5 eV for actinide elements and 8.0 cV 
for oxides. We used an fee structure for Pu, Am, and 
Cm with corresponding volume of each phases, while a- 
and /3- phases are used for U and Np, respectively. 

When the atomic interactions are turned off, the 
LDA-I-DMFT method reduces to the LDA method. Since 
there has not been a systematic study of how this method 
fares vis a vis the branching ratio, we include some LDA 
results in our study. In Fig. ^ we show the branching 
ratio of actinide elements obtained by the LDA method 
and compare to corresponding experimental data. We 
also show results for the three common angular momen- 
tum coupling schemes: LS (Hund's rule is dominant), 
jj (spin-orbit interaction is dominant), and intermediate 
coupling scheme (obtained from isolated atomic limit). 
The intermediate coupling scheme shows a good agree- 
ment with experimental results of late actinide elements. 
Am and Cm. In Ref. [Ill, the calculated value of Cm 
metal overestimates the experimental value. In our re- 
sult, it shows better agreement with experiment because 
we scaled the Slater integrals by 70%, which slightly 
moves the intermediate coupling scheme toward jj cou- 
pling scheme compared to the scaling of 80%. The con- 
ventional value 80% is reasonable approximation for very 
localized system such as oxides. However, smaller values 
are necessary for very itinerant systems. For example, 
for transition metal elements, the scaling factor is drasti- 
cally decreased for better description of spectroscopy^!. 
As pointed out by Refs. for a-U, and Ref. [l^ for a- 
Np, these metals deviate from the intermediate coupling 
scheme, which means these systems are delocalized and 
the effective spin-orbit strength decreases. 
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FIG. 3; The LDA+DMFT expectation value of the angular part of the spin-orbit interaction. The LDA+DMFT results are 
denoted by circles and rectangles for actinide elements and oxides, respectively. The experimental results obtained from EELS 
and XAS are denoted by thick dashes(l-8), because the number of electrons are not defined in the experiments. (1) a-U XAS^, 
(2) a-U EELS^ii^, (3) a-Np XAS^, (4) a-Pu XAS^, (5) a-Pu EELS^, (6) q-Pu EELS^, (7) 5-Pu EELSi^, (8) Am I EELSii, 
(9) Cm EELSii. 



The LDA results are in overall disagreement with ex- 
perimental data and the intermediate coupling scheme. 
From U to Pu the spin-orbit strength is much underesti- 
mated. Due to the overestimated band width of the LDA 
method, the spin-orbit strength of the actinide clement 
is considerably suppressed compared to the atomic cases. 
Only Am is well described by the LDA method due to 
its special configuration. There is an optimal spin-orbit 
stabilization of Am configuration, which gives clear 
splitting of occupied j = 5/2 and uoccupied j = 7/2 
states in LDA density of states. Our results show that 
the LDA description of the branching ratio is not proper 
for actinide elements, as shown in Fig. [21 

Figure [3] shows the branching ratio obtained by the 
LDA+DMFT calculation. Our results are in good agree- 
ment with experimental results of all actinide elements. 
The branching ratio of late actinides. Am and Cm, are in 
good agreement with the intermediate coupling scheme 
as experimental data. In well localized system, only sin- 
gle atomic multiplet configuration is occupied and the 
system can be well approximated by an atomic problem, 
which is the intermediate coupling scheme. On the other 
hand, as the system is dclocalized the / electrons be- 



come fluctuating between various atomic configurations 
and exchanging electrons with the surrounding medium. 
As a result, the spin-orbit strength is suppressed and the 
branching ratio moves from the intermediate scheme to- 
wards the LS coupling scheme. Note that the value for 
U and Np in Fig. [3] is further away from the intermediate 
coupling scheme than late actinides. Also note that this 
is not a competition between the spin-orbit iteraction and 
the Hund's rule coupling. This is the competition be- 
tween the spin-orbit interaction and the delocalization. 
If the system is fully dclocalized, the spin-orbit strength 
goes to zero and the branching ratio is just statistical 
value of 3/5. According to our plot, S-Pu also deviates 
from the intermediate coupling scheme, hence this sys- 
tem is partly itinerant. This deviation becomes larger 
for earlier actinides such as Np and U. 

Our results for U and Np slightly overestimate the 
spin-orbit strength measured in experiments. Here, we 
neglected the interference effect between different site in 
the description of core- valence transition as expalined in 
Eq. (4). In the impurity problem, we use the one crossing 
approximation method, which might overestimate the lo- 
calization of / electrons because it is based on the local 
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FIG. 4: Probability of the most occupied atomic ground state 
multiplet. The height of the peak corresponds to the fraction 
of the time the / electrons of the solid spends in the given 
atomic multiplet, denoted by the 5/ electron count Nf and 
the total spin J of the atom. 



atomic description. Also, the on-site Coulomb interac- 
tion U should be smaller than 4.5 eV for rather itiner- 
ant system, a-U and a-Np. By improving those effects, 
LDA-fDMFT can give more dclocalizcd branching ratio 
behavior. 

We also investigate the change of spin-orbit strength 
under pressure, which arc not available in atomic calcu- 
lation. The volume of a-Pu and Am IV is highly sup- 
pressed by 25% and 40% compared to the volume of 5- 
Pu and Am I, respectively. The difference of and the 
branching ratio between a-Pu and 6-Pu is almost negligi- 
ble and hard to detect, similar to EELS experiment a^^'^^ . 
On the other hand, there is a noticable change between 
Am I and Am IV. There is a transition from optimal spin- 
orbit stabilization for to optimal exchange interaction 
stabilization for configurationii. This transition in- 
duces a significant change in branching ratio between Am 
I and Am IV only with slight change of n^. 



The branching ratio shows systematic deviation from 
atomic calculation across the actinide series. Figure H] 
shows the probability of the most occupied ground state 
multiplet for each elements. The branching ratio of Cm 
and Am are very close to the atomic case, because the 
probability of the f (f) ground state is 95% (90%) in 
valence histogram^, which reflects almost atomic ground 
states. As the atomic number is decreased, the system 
becomes delocalized, and the spin-orbit strength is de- 
creased. In early actinides, the 5/ electrons are delocal- 
ized and they are distributed over various atomic configu- 
rations and not only in the ground state atomic multiplet. 
Note that 6-Pu already shows deviation from localized 
states, as discussed in the suppressed spin-orbit strength 
of (5-Pu. Under pressure, the 5/ electrons are delocalized, 
and the probability of atomic ground state multiplet is 
decreased as shown in Am IV and a-Pn. 

Finally wc also study actinide oxide systems PU2O3, 
Pu02, and UO2. These materials are very important 
because they are used in nuclear fuels for energy gen- 
eration. LDA-I-DMFT predicts that these materials are 
localized. As shown in Fig. [3] their branching ratios are 
consistent with the intermediate coupling scheme. The 
calculated n/ indicates that trivalent metal ion in PU2O3 
and tetravalent metal ion in Pu02 and UO2, which shows 
that these systems can be described by ionic system 
rather than metallic or covalent system. Recent experi- 
ments on the 4d 5/ transition in EELS show that the 
branching ratio of UO2 and Pu02 are similar to that of 
a-U and (5-Pu, respectively and it has been suggested 
that these actinide dioxides have covalent metal-oxide 
bonding, and the number of / electrons can become non- 
integer number—. However, unpublished EELS results 
of the 5d 5/ transition in UO2 and Pu02 suggest the 
bonding is indeed ionic with a near integer change in 5/ 
occupancy^^. So, further experimental work is required 
to investigate why the 5d and Ad spectra show ionic and 
covalent bonding nature, respectively. 
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